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Abstract 

> : 

■ We present results for the masses and decay constants of the light mesons 
in quenched QCD using the standard gluon action and a tadpole-improved 
SW-clover fermionic action to reduce discretisation errors. The calculation 
has been carried out at fixed volume and three lattice spacings corresponding 
to (3 = 5.7, 6.0 and 6.2. We make comparisons with the conventional SW- 
clover scheme. We use our results to extract continuum limits and to quantify 
the size of discretisation errors at smaller /3- values. 
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I. INTRODUCTION 



To extract reliable continuum results from lattice studies, it is essential to explore a 
range of /3- values. It is also valuable to employ a fermionic discretisation which has minimal 
discretisation errors. Here we investigate the tadpole-improved clover-SW action and we 
present preliminary results for the quenched spectrum at three lattice spacings. We compare 
our light hadron results at /3 — 6.2 with the conventional SW-clover approach and find no 
statistically significant evidence of further improvement at that /3-value. 

We use non-local sources to allow a more accurate extraction of ground state masses and 
couplings. Combining our results at all three /3-values, we can extract the continuum limit 
and assess the size of discretisation effects. We find that dimensionless ratios of observables 
have small discretisation corrections. We highlight the uncertainties remaining in extracting 
continuum decay constants because of the reliance on perturbative matching coefficients. 

II. THE SW CLOVER ACTION 

The SW fermion action is written in the form 

Sp = Sp (/?) KcS c \ oveT (1) 

with hopping parameter k. The quark fields q and q which appear in any observable must 
also be transformed (rotated). For an on-shell observable bilinear in the quark fields, the 
equation of motion can be used to express this as an improved observable given by 

0™ p = (1 + am q (l - z))q{T + zT®)q (2) 

where < z < 1. The size of the clover term in the action is determined by the coefficient 
c. At tree-level, c is equal to one. As demonstrated by Heatlie et al. 0, the terms of order 
a are removed for this choice of c. 

A mean field estimate of this coefficient was suggested by Lepage and Mackenzie 0. In 
this approach, the poor agreement of lattice perturbation theory with Monte Carlo simula- 
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tions of short-distance observables is explained by the presence of non-physical "tadpole" 
diagrams. One expects the power expansion for the gauge link 

U^x) = 1 + iagA^x) + 0(a 2 g 2 ) , (3) 

to be valid for sufficiently small a. However, tadpole diagrams introduce ultraviolet diver- 
gences which can partly cancel the factors of a 2 . A solution to this is to rescale the gauge 
link by a factor uq, such that U^{x) — * U^(x) = U^(x)/uo ■ The factor u is chosen so that a 
Monte Carlo calculation of a short distance observable using the rescaled gauge field agrees 
with the perturbative calculation. An appropriate choice is 

u = (ilH/ptaq) 174 • ( 4 ) 

Inserting this factor into the SW action changes the hopping parameter to k = kuq and 
the clover coefficient c from 1 to u^ 3 . A perturbative evaluation of c by Naik J1J gives a 
numerically similar value which is some support for the assumption of dominance of the 
tadpole diagrams. 

A. Computational Details 

Gauge configurations were calculated at (5 — 5.7, 6.0 and 6.2, using spatial volumes 
that are approximately the same. At each lattice spacing, propagators at 2 quark masses, 
straddling the strange quark mass, were calculated for each configuration — see table 1. 
Hadronic correlations were evaluated for each combination of such quarks, yielding three 
mesonic masses. Hence the extrapolation to the chiral limit could be explored using three 
hadronic masses. The configuration and propagators at (3 = 6.0 and 6.2 were calculated on 
a 320 node Cray T3D. The gauge configurations and propagators at (3 — 5.7 were calculated 
respectively on a 16 and 64 node i860 Meiko Computing Surface. For the transformation in 
eq. 0, we take z = which corresponds to no rotation in the propagators. 

The propagators were calculated using both non-local and local sources so that ground 
state and excited state contributions to correlators could be separated effectively. Local 
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sinks were employed. At (3 — 5.7, Jacobi smeared sources || were used, and correlations 
were measured with three combinations (neither, one, both) of nonlocal and local quark 
propagators in a meson. At (3 — 6.0, local and "fuzzed" sources were used |J. A fuzzed 
source can be written in the form 

3 

s^ed^) = £ F (y) s ( x -yi) > ( 5 ) 

i=l Vi=±r 

where F(y) is an SU(3) matrix representing the average of the fuzzed paths of gauge links 
from the origin to (y) (see fig. [I] for a schematic description). Fuzzing is less computer 
expensive than Jacobi smearing, as the only numerical work involves the averaging of the 
gauge links. It is gauge invariant and avoids time spent in gauge fixing. We use the fuzzing 
prescription with 5 iterations as in ref. || and we choose r = 6a at (3 = 6.0. An example of 
the advantage of fuzzing at the source is shown in fig. where an efficient extraction of the 
ground state vector meson is seen. 

III. FITTING 

We focus on the pseudoscalar and vector mesons. In order to evaluate accurately the 
central value and statistical errors, we fit to as many channels of data as possible. Further- 
more, to isolate the ground state for a wide range of time slices, we use a two-exponential 
fit. 

The two-point correlations can be expressed as a sum over transfer matrix eigenstates 
£<0|O 2 (f, t)Ot(0)|0) = £ cosh m 4 (t - N t /2) . (6) 

3 i m i 

Here c % = (0\O\H l ) where Oi^ are operators that have overlap with the ground and excited 
mesonic states \H°) and {H 1 ). These operators can be local or nonlocal, or have different 
gamma matrix structures. For both mesons, we fit simultaneously to all the available local 
and non-local source correlations of two-point functions. 

In the case of the pseudoscalar, we note that the fields P{x) = q(x)j 5 q(x) and A(x) = 
g(x)7°7 5 g(x) will create/annihilate the pseudoscalar, hence a simultaneous fit is performed 
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to two-point functions with the four combinations PP, AP, PA, AA. These contributions 
can be described in a factorising fit by the amplitudes (0\A\PS) and (0\P\PS), where \PS) 
is the pseudoscalar state. This provides tight constraints on the fit. 

We define the lattice decay constant for a pseudoscalar meson, fps as : 

(0\A L \PS) = m PS f PS , (7) 

where Al is the local axial current operator. The lattice / will be related to the continuum 
value / by including quark field rotations, setting the scale a and including the appropriate 
matching coefficient Za- One can determine fps directly from the factorising fit, without 
fitting to a ratio of channels as was used previously 0. An example of the result of such 
fits is shown in fig. |3|. 

Likewise, simultaneous fits for the vector meson can be performed using the operators 
V{ = q^q and V% = 97°7 l 9, and different smearing/fuzzing combinations. The lattice vector 
decay constant, fy is defined by 

£ £<V?(2, t)Vf(0)) -> -jf coshm y (t - N t /2) + A v coshm v (t - N t /2) , (8) 

j=l x JV 

where V{ is the local vector current. 

Despite the quite large statistics used in this calculation, the large number of different 
types of two-point functions requires that a full correlated fit must be carried out carefully. In 
particular, we employ the methods described by Michael and McKerrell || for the inversion 
of the correlation matrix. At /3 — 6.0 and 6.2, we model the correlation matrix with two 
exponents using the 5-diagonal approximation of ref. ||. At (3 = 5.7, we average all but 20 
of the largest eigenvalues of the correlation matrix. We fit to the largest range of t consistent 
with an acceptable x 2 - The errors quoted in this preliminary report are statistical only. 



IV. RESULTS: MASSES 
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A. Chiral Extrapolation 



The ground state pseudoscalar and vector masses are determined by the fits to two- 
particle correlations as described above. The issue of the rotations of the quark fields is 
irrelevant to mass determinations and will be discussed later. We define n cr - lt {j3) by the 
requirement m PS = for each lattice. At each lattice spacing, a linear extrapolation in 
the two squared pseudoscalar masses (where the quark masses are degenerate) was carried 
out. The results for k are given in table 2. A full fit using also the meson mass determined 
from the non-degenerate quark-mass case is possible, however, the fits gave statistically 
significant evidence for some curvature in mp s versus 1/ k. 

Without tadpole-improvement, there has been no agreement of non-perturbative lat- 
tice determinations of n CT it{f3) with one loop perturbation theory. For Wilson-like fermion 
actions, it is found that the tadpole-improvement factor of Uq makes k crit closer to the tree- 
level value of 1/8. For example, a determination of K CT \t using the c = 1 clover action |7j at 
(3 = 6.2 gave a result of n CT i t = 0.14315 ± 2 which corresponds to k crit = 0.12670 ± 2. 

For the tadpole-improved action explored here, we present in table 2 the comparison 
of the non-perturbative results for k cr \ t with the one loop perturbative calculation. Here 
we have chosen to use a coupling derived from a = — J^logSpiaq. Although the agreement 
is excellent at (3 = 6.2, the trend as (3 is decreased is not. This lack of agreement of K cr i t 
between non-perturbative determinations and one-loop perturbation theory at smaller (3 
will cause some uncertainty in normalisation when quark mass factors are included in the 
determination of the decay constants. 

B. m v — m PS 

The experimental value of the difference m v — m PS remains remarkably constant, at 
approximately 0.55 GeV 2 , over a wide range of masses. Quenched lattice simulations of rel- 
ativistic quarks have been in poor agreement with this result, and the discrepancy increases 
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as the quark mass is increased. In the heavy quark regime, such splitting has a leading 
order dependence on the chromo-magnetic coupling of the heavy quark with the gluon field. 
This is precisely the term which is amplified in the tadpole-improved clover action. While 
one cannot make a similar statement in the light quark regime, it is plausible that this 
modification of the clover term will also have a significant effect there. 

A comparison at (3 — 6.2 can be made with non tadpole-improved UKQCD data |7|||, 
based on 60 configurations. These configurations exist as a subset of the 184 configurations 
available now. A comparison of the results is shown in fig. ||. Unfortunately, there is no 
significant statistical difference between the two actions. We will soon have results for fuzzed 
sources in the tadpole-improved mass extraction which should lead to a reduction of errors. 

The behaviour at all three lattice spacings is shown in fig. [5| The scale (a" 1 ) is determined 
from the chirally extrapolated m p . The fall-off of the splitting with respect to m 2 PS is 
slightly faster than one would expect by comparison with the physical splittings m 2 — m 2 
and m 2 Kjf — m 2 K . However, until this difference is examined for quark masses around that of 
the charm quark, the advantage of tadpole improvement in this case remains unproven. 

C. Discretisation effects 

The dimensionless ratio of masses is known to have discretisation errors of order era for 
the SW-clover fermionic action. The tadpole- improved variant should preserve this result 
and is expected to improve on it by reducing the coefficient. We explore this by studying 
such ratios of observables against an observable proportional (to leading order) to a. 

Since the Wilson gauge action has inherent discretisation errors of order a 2 , it is advan- 
tageous to use this as a reference. We present ratios of m p to \f~K where K is the string 
tension determined by interpolating all existing data. In fig. |6|, we find that the ratio is 
constant (within statistical errors) at (3 = 6.0 and (3 = 6.2. There is however, a significant 
deviation from this constant for the ratio at (3 — 5.7. Since this behaviour is consistent with 
a 2 , it could be ascribed to the discretisation errors in either (or both) of the pure-gauge or 
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fermionic quantities. 

It is important to compare lattice results in the quenched approximation without chiral 
extrapolation if possible, since there are known to be anomalies in the chiral limit of a 
quenched theory. As suggested by Lacock and Michael 0, a dimensionless quantity can be 
constructed using this data without having to perform a chiral extrapolation: 

dmy , . 

J = mv d^- S (9) 

where the slope and my are to be evaluated at the vector mass satisfying my/mps = 1.8 
which corresponds to the K* meson. 

Our results are shown in fig. [7] and there is no evidence of any discretisation error since 
the values of J from each lattice spacing are equal within errors. Our result agrees with that 
from most existing quenched lattice data which are consistent || with a value of J of 0.37. 
This value is in disagreement with experimental data and, to substantiate this discrepancy, 
we will need a careful study of systematic errors, particularly at our largest /5-value. 



V. RESULTS: DECAY CONSTANTS 

As can be seen from eqs. [7| and [| the parameters fp$ and l//y can be easily obtained 
from the lattice data. However, to evaluate the physical quantities fpg and 1/ fv, one has to 
specify a suitable normalisation for the quark propagators and a renormalisation factor for 
the lattice operators. For fermion actions of the kind of eq. 1, several equivalent solutions 
have been proposed. To estimate the size of the systematic error induced by relying on 
perturbative estimates of the matching factors, we present here two interpretations of the 
action (eq. 1) with c = Uq 3 which result in relatively large differences. 

This systematic uncertainty is due to both our ignorance on the contribution of 0{a 2 s ) 
terms in the perturbative expansion, and to the possible different interpretations of the 
tadpole improvement procedure. 
Normalisation I 
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In this case we interpret the clover coefficient c = Uq 3 as a consequence of the tadpole 
improvement, which implies a redefined hopping parameter (k) and a gauge field rescaling. 
Because — > U^/uq, the covariant derivative (or the lattice bare mass for the unrotated 
case 2 = which we are using) which appears in eq. 2 has to be divided by uq. We 
choose, however, to adopt a slightly different normalisation, which amounts to using the 
renormalized quark mass, 

am q = I f i - 1) , (10) 



2 \ Kj k>c 

instead of the bare mass since the "rotation" factor 1 + am q is then relatively closer to unity 
and possibly better reproduced by perturbation theory. This is possible because formally 
the difference between the bare mass and the renormalised mass only contributes at order 
aa which does not effect the improvement. 

The renormalisation constants to be used in this case differ from those corresponding to 
the SW action, and have been computed in ref. [fTOH , and are 



Z A = 1 - 0Ala s + 0{a 2 s ) , 

Z v = 1 -0.58a s + O(a 2 s ) , (11) 

where we choose to use a s = — ^ log Spiaq as an appropriate improved estimate of the 
coupling. There is an inherent systematic error from the choice of scale for this coupling - 
this can only be resolved by a two-loop calculation. 
Normalisation II 

Let us consider the results one obtains using the action of eq. 1 with c = Uq 3 but without 
attributing to such a coefficient the meaning of tadpole subtraction. Because of that both 
the hopping parameter and the normalisation of the propagator itself keep their original 
values. Moreover, since the perturbative expansion of the clover coefficient is 1, up to terms 
of 0(a), the renormalisation constants are the same as those calculated for c = 1, given in 
ref. 0. In this case we used the "unrotated" operators corresponding to taking z = in 
eq. 2, where m q is now interpreted as the "bare" lattice quark mass. 
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In order to compare with other work, at each /3, we extrapolate the lattice decay constants 
to the chiral limit assuming a linear behaviour with m 2 PS . The dimensionless quantities 
f w /y/K and are shown for each normalisation in figs. ^. 

The results for f n / \J~K for both normalisation schemes are consistent with having a 
common continuum limit which agrees with the experimental value. As (3 is decreased, the 
difference between normalisations I and II increases. This largely arises from the fact that 
the perturbative estimate of m q is not in agreement with the numerical value (it is therefore 
due to large unknown contributions at 0(a 2 )), and this causes a poor renormalisation in 
the case II. An insight into this problem could be achieved by calculating the "rotated" 
operators of eq. 2 with z — 1, which would permit us to eliminate the uncertainty in the 
renormalisation of the quark mass. Such a study is currently in progress. Ideally, however, 



one should measure the matching constants non-perturbatively, as outlined in refs. [ I I and 



0- 

The c = 1 results || for J^/yK at (3 — 6.2 are seen to lie significantly lower in fig. |8|. 
A careful study of systematic errors is needed to explore this discrepancy further. Since we 
have used different rotation factors (z = for c = Uq 3 and z = 1 for c = 1), the effect 
of this difference on the discrepancy can also be determined by a study of z = 1 for the 
tadpole-improved case. 

Fig. |9| shows for l// p a rather similar picture of agreement between the continuum limit 
of the tadpole-improved schemes and experiment. In this case the c = 1 result is also in 
agreement. 



VI. CONCLUSIONS 

The goal of lattice investigations of QCD is to find a formalism which has small discreti- 
sation errors. This has to be studied by obtaining results for a wide range of observables for 
several lattice spacings. In this paper we have outlined preliminary results on the spectrum 
and couplings from a tadpole-improved SW clover action. We find that, within our present 
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statistical errors, the results at (3 — 6.0 and 6.2 are consistent while those at (3 — 5.7 show 
signs of discretisation effects. This is encouraging since it suggests that meaningful results 
can be obtained at 6.0 which corresponds to a -1 > 1.8 GeV. To push the frontier down to 
5.7, it seems plausible that an improved gauge action will also be needed since the ratio of 
m(0 ++ )/\/K evaluated using the Wilson gauge action is known to deviate significantly at 
j3 = 5.7 so that there are significant order a 2 discretisation errors here. 

We also emphasize that continuum extraction of decay constants and matrix elements 
is beset with ambiguity if a perturbative matching is used. This arises in part because 
the critical hopping parameter is not accurately given by 1 loop perturbation theory. A 
non-perturbative determination of these matching coefficients is essential for fully reliable 
continuum predictions. 
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TABLES 






Ng x N t 


Number of 
configurations 


K 


Smearing 


Clover 
coefficient 


5.7 


12 3 x 24 


480 


0.13843, 0.14077 


Jacobi & Local 


1.5678 


6.0 


16 3 x 48 


130 


0.1370, 0.1381 


Fuzzed &: Local 


1.4785 


6.2 


24 3 x 48 


184 


0.1364, 0.1371 


Local 


1.4424 



TABLE I. Input parameters for simulation. 



p 


^crit 


~ lLoop 
K crit 


am p (K CTit ) 


a-), [GeV] 


5.7 


0.12344 +2 

- 2 


0.1214 


0.678 +7 

- 7 


1.14 +1 

- l 


6.0 


0.12218 +3 

- 2 


0.1219 


0.405 +13 

-12 


1.90 +6 

- 6 


6.2 


0.12206 +1 

- 2 


0.1221 


+12 

0.302 

-14 


2.55 +12 

-10 



TABLE II. Output results of simulation. Non-perturbatively derived values of R cr \t compared 



with a 1-loop calculation using tadpole-improved perturbation theory. The vector masses extrap- 
olated to the chiral limit (m^ = 0) and the resulting scale, taking m p = 0.7699 GeV. 
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FIG. 1. Schematic guide to fuzzing. The source is composed of six points around the origin. 



The dotted lines indicate typical paths that will be used in the average. 
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FIG. 2. The effective mass of the vector meson versus time separation t. The data are at 
(3 = 6.0 and k = 0.1381, using local (diamonds) and fuzzed (crosses) sources. The fits shown are 
two state expressions. 
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FIG. 3. Effective mass plot from the local and smeared combinations of the operators A and 
P at = 5.7 and k = 0.13843. 
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FIG. 4. The comparison of the mass combination (amy) 2 — (amp) 2 at [3 = 6.2 from SW-clover 
and tadpole-improved SW-clover fermionic actions. 



14 



0.8 



0.6 



> 
a; 
O 



0.4 — 



0.2 



0.0 



o /S = 6.2 
o ,3 = 6.0 
□ /S = 5.7 



0.0 



0.5 



1.0 



FIG. 5. The comparison of the mass combination (my) 2 — (mp) 2 between lattice results and 
experiment. 
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FIG. 6. The dimensionless ratio of the p-mass to the square root of the string tension. The 
"experimental" point at a = corresponds to (770MeV/440MeV). 
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FIG. 7. The dimensionless quantity J = m v dmy / drn? PS from the different lattice /?- values 
compared to the experimental value. 
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FIG. 8. The dimensionless ratio of f w to the square root of the string tension. The "exper- 
imental" point at a = corresponds to (132MeV/440MeV). The points are displaced slightly for 
clarity. 
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FIG. 9. The p meson decay constant 1/ f p versus the lattice spacing a. The points are displaced 
slightly for clarity. 
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